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Mixing in open incompressible flows is studied in a model problem with inhomoge- 
neous passive scalar injection on an inlet boundary. As a measure of the efficiency of 
stirring, the bulk scalar concentration variance is bounded and the bound is shown 
to be sharp at low Peclet number. Although no specific flow saturating the bound at 
high Peclet number is produced here, the estimate is conjectured to be approached 
for flows possessing sufficiently sustained chaotic regions. 



I. INTRODUCTION 

Since the pioneering work of L. N. Howard [1] , applied physicists and mathematicians have 
found that rigorous bounds are a straightforward way of obtaining useful estimates on fluxes 
of quantities transported by complex incompressible flows, from laminar to fully developed 
turbulent flows. In recent years the co-authors and collaborators [2-7] and others [8-12] used 
similar ideas to bound mixing efficiencies for processes described by the advection-diffusion 
equation with inhomogeneous passive scalar sources and sinks. One major insight emerging 
from these studies is the distinction between transient mixing problems, where concentration 
inhomogeneities are provided by initial conditions, and (statistically) steady mixing in the 
presence of sustained sources and sinks. In the former case the flow's shearing and straining 
properties, i.e., the stretching and folding of material lines, serve to enhance the effect of 
molecular diffusion and increase the rate of scalar concentration variance dissipation. On 
the other hand, in the latter case the spatial variation of the concentration is best suppressed 
by direct transport of material from sources to sinks and vice versa. Then the maximum 
possible enhancement of molecular diffusion by stirring is Hmited not by small scale features 
in the flow, but rather by small scale properties of the sources and sinks [3-5] . 

In this paper we extend these techniques to the case of fluxes at the boundaries rather 
than body sources as previously described. This is not just a mathematical detail: the 
problem changes significantly. Notably, the scalar inhomogeneities are transported away 
from both source and sink regions so the transient mixing efficiency of downstream regions of 
the stirring flow are crucial for efficient mixing. At the same time the source-sink structure 
defined by the distribution of fluxes on the boundary plays a key role in the maximum 
possible mixing efficiency. In the remainder of this paper we define a specific scenario where 
these issues can be addressed and study it via general analysis and particular examples. 
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FIG. 1: Schematic of the flow domain F. The scalar flux is imposed at x = and x = h, and there 
is no flux on the sidewaUs. 

II. A BOUND ON THE VARIANCE 

The suppression of concentration variance is a reasonable measure of the efficiency of 
mixing since in the absence of sources it decays due to diffusive mixing and it vanishes when 
complete homogenization is achieved. The variance is also a useful measure in the context 
of open flows [13-15]. We use techniques developed in [2-6] to bound the time-averaged 
variance of the concentration of a passive scalar. We shall use both the variance and its 
square root, the standard deviation, but final results will usually be quoted in terms of the 
variance. 



A. The model and definitions 

Consider a fixed domain F, in which flows a fluid with incompressible velocity fleld u(x, t). 
The fluid transports a diffusing passive scalar with concentration 6{x,t) which is stirred to 
some degree of homogeneity within F. The concentration 6{x, t) in the domain F obeys the 
advection-diffusion equation 

dte + u-we^ KAe. (i) 

The flux of 9{x, t) is deflned by 

F:^ue-KWe. (2) 
At the boundary dV of the domain, we specify the normal flux F ■ h as 

{F ■h){0,y,z) = J^Q{y,z), {F ■ h){h,y, z) = J^hiv, z), • n)side walls = 0, (3) 

where the sidewalls are the curved sections in Fig. 1. We denote the horizontal flux Tix) := 
x-F{x), and we will often write J^xiv, z) to mean J-'(a;, z)^ as in the boundary condition (3) 
above. Note that there is considerable literature, and a long running debate, regarding the 
appropriate boundary conditions to use for open-flow systems [15-19]; we use a prescribed 
normal flux for mathematical expediency. 

Multiply the advection-diffusion equation (1) by an arbitrary smooth test function <^(a;, t) 
and integrate: 

J ^dtOdV + J ■ {ue)dV ^ K J ^AedV (4) 
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where dV^ = dxdy in 2D and dxdydz in 3D. To obtain (4), wc have also used the incom- 
pressibihty condition V • w = 0. Next foUow a few integration by parts, 



(5) 



dtj OipdV- J 9dtipdV + J V-{uip9) - j OV-{uLp) dV 

= ^(^J V ■{^VO)dV - J V ■{9V^)dV + J 9A^dV^ . 

We collect on the right the boundary terms, 

dtj ecpdV- J e{dtcp + V ■{ucp) + KAcp) dV ^- J V ■{cpF + K0Vip)dV , (6) 

and turn them into a surface integral, 

RHS = -/" {ipF + k9V(p) -nda, (7) 
Jar 

where da = dy in 2D and dydz in 3D. On the sidewalls, we have F • h — because of 
boundary conditions (3), so the flux term in (7) vanishes. Since we haven't solved the full 
problem (1), we don't know 9 on the boundary. To remove the direct 9 dependence in (7), 
we restrict attention to test functions ip with 

Wip-h^O on dr. (8) 

All that remains of the integral (7) are the surfaces at a; = and x = h, and Eq. (6) is then 

dt [ 9ipdV + {9C*ip) = [ ifJ^oda- [ ^J^^da, (9) 

with da = dy in 2D and dyd^; in 3D, and dVj. is the cross-section at constant x. We also 
defined the formal adjoint of the advection-diffusion operator, 

jr*cp ■= -dt(p - V • (w - kA(p . (10) 

Equation (9) is the fundamental 'integrated equation' that we will be using subsequently. 
Setting = 1 in (9), we get an equation for the conservation of total 9, 



dt 9dV= Joda- / J), da (11) 

which of course says that the change in J^9dV is the net difference in fluxes at x = 
and X — h. For the remainder of the paper, we shall consider the case where 



Toda= / J'/.da^O (12) 



so that the total amount of ^ in F is always 0. 
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B. Lower bound on the vsiriance 

We define angle brackets and overbars as 



if) lim U dt [ f{x,t) dV, f := lim U dt [ f{x,t)da. (13) 

(The limits are assumed to converge.) After time-integrating Eq. (9) the term {dt{(fi9)) 
vanishes and we have 



{eCif)^ {ipF + ned^ip) - (ipT+Ked^if) . (14) 







The notation /|^ signifies / evaluated at x = c. The conservation equation (11) im- 
plies {9) = 0, so the spatially-integrated variance of the concentration field is (6'^). We 
shall omit the qualifier 'spatially-integrated' most of the time. 

Using the Cauchy-Schwartz inequality, we can extract from (14) a bound on the standard 
deviation, 

- ■ ^^^^ 

The right-hand side of the bound (15) contains only known quantities. The bound says 
that for given functions u, J-q := J-'Iq, J-'h '■= J^l^, and diffusivity k, the passive scalar must 
exhibit a minimum level of fluctuations. Since better mixing requires decreasing fluctuations, 
Eq. (15) is a rigorous bound on the potential efficiency of a mixer. In practice, one tries to 
find a function (/? to make the bound as sharp as possible, that is to maximize the right-hand 
side of (15). The optimal (p will, in principle, depend on all the problem parameters. 

Observe that adding a constant c to </? changes the numerator in (15) by a factor c {J^o — 
Jti), which vanishes because the total amount of 9 within the domain is conserved (after the 
long-time averaging). Hence, adding a constant to (p does not affect the estimate. 

If we have a semi-infinite 'channel,' with h — )■ oo, then the scalar field at h will become 
completely homogeneous due to diffusion (6* — )■ 0) and its flux will necessarily vanish. We 
then have 



Equation (16) is ideal for this channel configuration, since we only require knowledge of the 
flux at a; = 0. (Alternatively, we could also obtain (16) by specifying = aX x = h.) We 
shall focus exclusively on the semi-inflnite channel conflguration [h — )■ oo) for the remainder 
of the paper. 



C. A bound for all flows 

The bound (16) requires complete knowledge of the velocity fleld u. However, we can 
turn it into a 'global' bound over all allowable velocity flelds. For simplicity, from now on 
we restrict attention to a uniform, rectangular cross section of dimensions Ly and L^, and 
to a time-independent velocity fleld. 

First we bound the numerator in (16). For a time-independent test function (f, 

(|£V|2) = {\u ■ Vcp + kAcpI'') < 2{\u ■ V(^|') + 2«2(|A(^|2). (17) 
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We bound the first term in (17), 

2(|wV(^H<2/ |u|2|V(^|2dx 

JO 

POD 

< 2 / ^suplVv^l^da; 
Jo y,z 

POO 

<2sup\u^ sup \\/ ip\^ dx ^2U^L'^-^Ci 
X Jo y,z 

where we defined the mixing fiow velocity scale C/ > by 



:= L^-" sup |wP (18) 

X 

and the length scale L :— ^jLyL^ in 3D, or L := Ly in 2D. We also defined the dimensionless 
coefficient 

POO 

CiM:=L / sup|V(/?|Mx. (19) 

Jo y,z 

We define another dimensionless coefficient 

c,W} :^ L^-'^ {\A<p\') (20) 
which, using (17), allows us to write 

(|£Vr) < 2U^L''-^ ci + 2«:2l'^-4 = 2^1''-^ (ci Pe' + C2) (21) 
where the Peclet number is defined as 

Pe UL/k. (22) 

Now we examine the numerator in (16). The fiux has units of velocity (consider- 
ing without loss of generality that d is dimensionless), so we can define the dimensionless 
coefficient 

c^W.Fo}:^U-'L'-' . (23) 

Combining the results we can write the bound as 

^ ' - V2(ci Pe^ + c) ^ ^ 

where now the properties of the stirring flow only appear as the scale U in Pe. For small 
or large Pe, we can improve the bound a little by neglecting the appropriate term in (17), 
which allows us to leave out the factor of 2 in that equation. We then get the asymptotic 
estimates 



..2^1/2 > \L^'\c,/^^)Y>e, Pe« 1; 

^\L'^'\c,lJcll Pe»l. ^ ^ 



For convenience of display, we combine the asymptotic bounds (25) into the same form 
as (24), 

^g2^i/2>^d/2 C3Pe ^ Pe«lorPe»l, (26) 
V ci Pe^ + C2 
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where it is understood that this estimate is only vahd asymptotically for small or large Pe. 
For our purposes, the error made at intermediate Pe will be small. 

The small Pe bound says that the standard deviation can become small in that case: 
this is the diffusion-dominated limit. For large Pe, the lower bound saturates and becomes 
independent of Pe: we shall examine this behavior more closely in Section IV, and see if it 
can be realized in practice. 



D. Choosing the test function ip 

To get usable numbers out of (26) we need to choose ip for a given J^q. We specialize to 
two dimensions {d = 2) and 

My) ^ UH sin qy, (27) 

where the constant is a dimensionless measure of the flux. We consider periodic boundary 
conditions in y, which for our purposes is equivalent to the no-flux conditions on the sidewalls 
in (3). The wavenumber q = Q/L, where Q is an integer multiple of 27r and L = Ly. 

To get the sharpest bound, our goal is to maximize the right-hand side of (26) as much 
as possible. Hence we need a test function (p with a large projection on Tq, to maximize C3. 
The most expedient approach is to take a separable form for (f with the y dependence 
mimicking J-q^s: 

if = $(a;) singy, $(0) = 1. (28) 

We can then easily compute 

C3 = U-^L-^ $(0) UH I sin^ qy dy = \H. (29) 

Jo 

and 

/■oo 

C2 = L2(|A(^|2) = L2(|($"-g2$)singy|2) = IlM - q^^'' dx . (30) 

The constant ci is more problematic, since we need to evaluate the supremum inside the 
integral: 

/■oo /■oo 

ci = L / sup |V<^|^dx = L / sup (($')^sin^g|/ g^$^cos^gy) dx. (31) 
Joy Joy 

To find the supremum, take the derivative and set to zero, 

(($')' sin^ qy + cos" qy) = {i^'Y - ?'^') Q sm 2qy = 0. (32) 

ay 

We need only consider the extrema at qy — and 7r/2, since sm2qy function is 7r-periodic, 
which leads to 

/■oo 

ci = L I max (5^$^ , $'2) da; . (33) 
Jo 

Now we make specific choices for $. Consider first the limit of small Pe. In that case we 
ignore Ci and try to minimize C2. The Euler-Lagrange equation for (30) is 

- 2g2$" + q^^ = 0. (34) 
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The appropriately integrable solution that also satisfies $(0) = 1, $'(0) = is 

= (1 + qx) e"^^ (35) 

which gives C2 = {qL)^ = Q^- 

Conversely, for large Pe we ignore C2 and focus on Ci. In that case, we might like to use 
the function $ = e"^^, for which = — so the two choices for the max in (33) are the 
same, and Ci = ^{qL) = \Q for the constant. The only problem is that this choice for $ 
does not satisfy the boundary condition $'(0) = 0. This is easily remedied by inserting a 
small boundary layer near a; = 0, and letting the size of the layer go to zero. For large Pe, 
the correction incurred can be neglected. 

We can combine our two asymptotic results as in (26) to conclude 

7^(g')^/'> , Pe«lorPe»l. (36) 

2y^|gpe2 + g3 

This estimate is rigorous for large and small Pe, and we expect that error made at inter- 
mediate Pe is typically small. This is illustrated by the example developed in the next 
section. 



III. EXAMPLE: THE BARBERPOLE FLOW MIXER 

As an illustration of typical features of open-fiow systems, we study a toy problem that 
can be solved exactly. We will then compare the exact standard deviation to the bound (36). 
We consider the 2D spatially uniform flow 

u — ux + vy , (37) 

where u and v are constants. We we solve the steady advection-diffusion equation for 9{x), 
where x = {x,y) & F — [0, oo) x [0, Ly] , with periodic boundary conditions in y and the 
flux (2) specifled as 

x-F(0,y) = J'o(?/), F{x,y)^0 as x ^ oo. (38) 

For the case of pure advection, with /t = in Eq. (1), the solution is 

9{x,y) ^ u~^To{y - xv/u) . (39) 

The concentration has a 'barberpole' pattern (see Fig. 2): the imposed distribution at a: = 
is twisted around the periodic direction y, with the tilt angle 7 determined by the ratio of v 
and u. The cross-sectional variance 0"^ /Ly is independent of x, and there is no mixing, as 
expected. 

A. Advection and diffusion 
Now we include diffusion in Eq. (1). We flrst Fourier-expand in y, 

e{x, y) = Y^ Oqix) + c.c. (40) 

9 
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FIG. 2: Concentration field 9{x) for the barberpole flow (parameters values are H = 1, Pe = 100, 
7 = 37r/8, L = Ly = 1). 



where q is an integer multiple of 27i/Ly. In the steady case, the advection-diffusion equa- 
tion (1) for the amplitude of the vertical mode q becomes 



kO" -uO' - q{Kq + iv) 



0. 



(41) 



Using the length scale L = Ly and velocity scale U = y/u"^ + v"^ (consistent with the def- 
initions in Section II C), we define the rescaled variables X = x/L, Y = y/L, Q = qL, 
and 0(-^) = 6q{LX) to get the dimensionless ODE 



e"-Pe cos7e'-Q(Q + iPe sin 7)6 = 



(42) 



where the pitch angle 7 is defined by tan 7 = v/u and the Peclet number Pe is as in Eq. (22). 
This has solution 

e(X) = 6(0) e-('^+^'^)^, X > 0, (43) 
where we defined the spatial decay rate by 



u = a — ^Pe cos 7 



(44) 



with 



«= VH^' + (iPe'cos2 7 + Q2)) 



i(i?2_(lPe2cos2 7 + Q2)), 



i?2 = ^g2pe2sin2^^ QPe2cos2 7 + g2)^. 

We have u > since a > (| Pe^ cos^ 7 + Q'^) > ^Pe cos 7; hence, the solution (43) decays 
with X, as required. Using the form (27) for the flux J-'o{y) and the boundary condition (38), 



K 



Toiy) = UH sin qy = T{0, y) = - (Pe cos 7 + (z/ + iw)) 9(0) e'"^ + c.c. 



(45) 



which we can use to solve for 



9(0) 



1 



HPe 



2i Pe cos 7 + 1/ + icj 



(46) 



Figure 2 shows the concentration field pattern for the parameter values H = 1, Q = 27r, 
Pe = 100, 7 = 37r/8, L = Ly = 1. The tilted stripes make an angle 7 with the horizontal, 
and the intensity decays in space at a rate u. 
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FIG. 3: Decay rate v from Eq. (44) for the barberpole flow, as a function of Pe. For any value 
of 7 < 7r/2, the decay rate ultimately goes to zero like Pe~^ for Pe ^ 1. However, for 7 = 7r/2 the 
decay rate grows as Pe^^^. 

For small Pe, we have the expansions 

i/ = Q-lPecos7 + 0(Pe2), w = ^Pe sin7 + O(Pe^) . (47) 

To leading order in Pe, the spatial decay rate v is independent of Pe and the pitch angle 7. 
This is because in the absence of flow (or for a weak flow) the decay scale length is set by g, 
and the diffusivity completely drops out. 
For large Pe, we have the expansions 

i^ = Q2sec^7Pe-^ + 0(Pe-3), w = gtan7 + 0(Pe~2) . (48) 

These are valid for 7 not too close to 7r/2. Figure 3 shows the spatial decay rate v as 
function of Pe, for several values of 7. For large Pe, the decay rate ultimately decays 
as Pe~^. However, for 7 = the decay rate grows as Pe^/^ This case corresponds to no 
horizontal flow: the vertical flow is very effective at pushing the source onto the sink, leading 
to fast spatial decay. But the flux into the mixing region is then entirely due to diffusion, 
so this is not a very practical situation. We will discuss this case further in Section IIIB. 
Another striking feature of the curve in Fig. 3 is that some curves have a local maximum. 

The enhancement to the decay rate is related to shear or Taylor dispersion, in the sense 
that it arises from the tilting of concentration contours. However, here it is not due to shear, 
but to the cross-velocity v. 

B. Computing the variance 

For a single Fourier mode, the variance integrated over a cross-section is 

W= 9^dy = 2L\e{X)f = 2L\e{0)fe-^'''', (49) 
Jo 
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FIG. 4: Comparison of the standard deviation of the concentration field for the barberpole mixer 
(sohd lines) and the global (lower) bomid (26) for all flows. For small Pe, all the curves are linear 
in Pe. For large Pe, the barberpole solutions scale as Pe^'^^, unless 7 = 7r/2 for which the scaling 
is PeV4. 



and the total variance in the domain F = [0, 00) x [0, Ly\ is 

/•CO T 

{d^)= 2L|e(0)|'e-2'^^da; = L2-|e(0)|', (50) 
Jo ^ 



or in term of the standard deviation 



-^{9^)^'= ^ ^ (51) 
-t^-n 2^^^ (Pe C0S7 + z/)2 + 



For small Peclet number, the standard deviation is 



_^//)2\i/2 Pe Pe2cos7 



,^.>"=^^-^ + 0(Pe3). Pe«l. (52) 
whilst for large Pe it has the form 

-^(^^2)i/2 = ^v/p^^ + 0(Pe-3/2), Pe»l. (53) 

If 7 = 7r/2, then the O^Pe^''^^ term in the large Pe expression diverges, and we have 

Z^(^V''^(^)"%0(Pe-/.). Pe (54) 

These scalings are evident in Fig. 4, which compares the exact standard deviation (51) to 
the bound (36). For small Pe the exact standard deviation and the bound are both linear in 
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Pe. For large Pc the lower bound flattens out and becomes independent of Pe. The closest 
the exact standard deviation for the barberpole flow comes to the theoretical lower bound 
is for 7 = 7r/2, since then the exact result scales like Pe^''^. For any other value of 7 < 7r/2 
the Pe^^^ scahng holds for a short while, but eventually rises as Pe^^^ when Pe cos 7 becomes 
large. Hence, the 'optimal' barberpole flow is one for which the horizontal flow u vanishes 
(7 = 7r/2): in that case all of the velocity field is devoted to pushing the source onto the 
sink, in a manner similar to the optimal flows found for closed systems [4, 7]. However, 
unlike the closed systems this flow does not scale optimally: the bound suggests that there 
could exist a flow producing a shallower scaling. We will discuss the possibility that such 
flows exist in Section IV. Note also that as mentioned in Section HI A the flow with 7 = 7r/2 
is very singular: the scalar flux into the domain is entirely due to diffusion, so there is no 
mass flux into the mixing region. This is not a very realistic concept on which to build an 
industrial mixer. 



IV. EXAMPLE: A MEANDERING FLOW MIXER 

The barberpole flow of Section HI is a nice model because of its exact solvability but 
it has some undesirable pathologies. Most importantly it is hard to relate it to practical 
situations where a channel flow enters and then exits a well-deflned 'mixing region,' where 
mechanical stirring occurs and most of the mixing takes place. A better example would 
involve a flow that varies in x. 

As a more general version of the constant barberpole flow (37), we now consider the case 
for which the y component is a functions of x, 

u{x) = ux + v{x) y (55) 

with u > 0. The domain and boundary conditions are as in Section HI. This is a 'mean- 
dering' velocity fleld that changes its direction with x but never backtracks. Once again 
we Fourier expand as in (40) and focus on a single mode 9q{x). The equation to solve is 
still (41) but with v{x) nonconstant. Unlike the constant flow case, we cannot solve the 
equation in full generality, but since we are mostly interested in the small k, case the WKBJ 
method allows us to get an explicit approximate solution. We assume the standard WKBJ 
ansatz 

eg{x) = ^,(0) exp(>S(x)) = OgiO) exp{So{x) + kSi{x) + ---) (56) 

where S{0) — 0, and to save notation we left out a q subscript on S. After inserting into (41) 
and collecting terms in k, we can easily solve for the flrst two terms, 

So{x) = -i^ rviOdL (57a) 
^ Jo 

^i(^) = -— -4K^)-^(0))-4 ^^'(0d^ (57b) 

We will not be flnding higher-order terms, and thus decree S = Sq + Si, it will prove more 
convenient to decompose this into real and imaginary parts as S{x) — S,.{x) + iSi{x), with 

S,{x) = (^u'x + J\\OdC^ , (58a) 

S,{x) = -^ rv{^)dC-^{v{x)-v{0)). (58b) 
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FIG. 5: The domain for the meandering flow: within the mixing region the flow takes the form (55), 
and for x > Lmix it has the simpler, constant form (u, 0) . 



Now we need to apply boundary conditions at a; = 0, where the flux is 

^■(0, y) = ^^(0) {u - kS'{0)) e^'^y + c.c. = UH sin qy (59) 
where we used S{0) — 0. We use this to solve for 

^''°) = 2i(„!fg-(0)) - 

Note that the derivative in S' simply annuls the integrals in the definition of 5", so this 
pref actor only involves the unintegrated velocity field at a; = 0. The variance integrated 
over a cross-section for a single Fourier mode is 

W= 9^dy = 2Ly\9g{Q)fe^'^^^''\ (61) 
Jo 

with 

1^,(0)1 = . (62) 

Thus S^{x) characterizes the spatial decay of the concentration field with x. The total, 
spatially-integrated variance is 

/•oo 

{9^)^2Ly\eg{0)f e2^--(^)dx. (63) 
Jo 

To make further progress we will now assume that we have a mixing region localized in x 
of length Lmix, extending from < x < L^i^ (see Fig. (5)). Before and after that region we 
have V — as well as v'{0) — 0. In that case S'{0) = 5*^(0) = —nq^/u. As for the integral 
in (63), it splits as 

/ e2^^(^)dx=/ e2^^(^)dx+ / e-^'^'^^/^dx. (64) 

Jo Jo Jl^Ix 

Since ^ n is assumed small, the first integral is roughly equal to Lmix and the second 
integral can be done exactly, yielding 



^0 



2^5^ 
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We already assumed that Kq^L^ii^/u was small in approximating the first integral by Lmix, 
so we do it again in the exponential term to find after expanding 







e25.(.) ~ At . (66) 



Neglecting the k, in the denominator of (62), we combine it with (66) in (63) to find 



Recall the definition of U, Eq. (18): outside the mixing region we have \u\'^/Ly — v?^ but 
inside we have \u\^ jhy > because of the extra stirring motion. Hence, we put u — cj^ C/, 
with C4 > 1, and obtain 

^(9y/^.i|PeV^ (68) 

where the length scale L = Ly and Q = qLy. For a constant fiow, we have C4 = 1 and we 
recover the asymptotic form (53) with 7 = 0. The constant C4 ^ 1 measures the vigor of 
stirring within the mixing region as compared with the mean fiow. 

The standard deviation paradoxically appears to increase with C4: this is an artifact of 
our dimensionless scaling, which was natural for the bounding approach of Section II C but 
is less appropriate here. If instead we normalize the fiux as J-'o{y) = uHoo sinqy, where 
now measures the fiux normalized by the infiow/outfiow velocity u, we find 

(^2)^2L,^-^ = L,^, (69) 

or 

Pe.:^«^. (70) 

So if we use a Peclet number based on the infiow, the variance is completely independent 
of the stirring fiow, at leading order in Pcoo- The stirring is too weak to change not just 
the Pe^^ scaling, but also the prefactor. This might seem surprising but is actually rather 
obvious: we assumed in calculating Eqs. (68) and (70) that the small diffusion dominated 
even in the mixing region. This is justified since the fiow is 2D and steady, so chaotic 
advection cannot occur [20]. Hence the finite mixing region does not decrease the variance 
enough, and the total variance is dominated by the infinite region x > Lmix- 

The scaling (68) for the general 'meandering flow' shows that there is no hope of ap- 
proaching the lower bound Pe" scaling with such a flow. Let us assume that the fluid 
gets well-mixed in the mixing region, in a diffusion- independent way as occurs for chaotic 
and turbulent mixing. Then the first term in integral (64) yields a diffusion-independent 
constant Lmix, and the second term can be neglected. The variance Eq. (67) becomes 

7/2 //2 

{e')c^2Ly—^L^^, (71) 

which leads to the standard deviation 

^(^^)V2:.1(L^,,/L)V2. (72) 
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This indeed exhibits a Pe-independent scahng. We emphasize that the difference be- 
tween (72) and (70) is that the former assumes that the length of mixing region required to 
reduce the variance to a certain level becomes independent of the diffusivity as k, becomes 
smaller, which is a standard assumption of turbulent and chaotic mixing. (In fact, typically 
there will be a weak logarithmic correction to Lmix, having to do with the length needed to 
create small scales.) In any case it remains an open challenge to produce an explicit flow 
that realizes such optimal mixing. 

V. SUMMARY 

We have proposed a framework for the theoretical investigation of questions of mixing 
efficiency in open flows with passive scalar tracer concentration variations sustained by 
sources and sinks at an inlet. The measure we used to gauge the effectiveness of the stirring is 
the spatio-temporal variance of the scalar concentration, a quantity that can be analyzed via 
the advection-diffusion equation. We derived lower bounds on the variance that are observed 
to be sharp at low Peclet number (weak stirring) and to capture the qualitative trend at 
high Peclet number (strong stirring) for a particular model problem, the "barberpole" ffow. 
Concentration variations introduced at an inlet and subsequently swept downstream are 
likely most effectively stirred (for the mixing criterion adopted here) by ffows that, in a 
downstream moving frame, provide the kind of chaotic mixing known to work well for 
transient mixing problems. It remains an open question to determine if such a ffow can 
be designed or characterized in detail even for the simple idealized model considered here. 

Many other questions remain open as well. How does the mixing efficiency of a given 
flow vary as the inlet source varies? How does the mixing efficiency depend on the length 
scales we choose to focus on? For example the gradient variance or some other norm might 
just as well be invoked as the gauge to measure mixing efficiency. This issue is known to be 
important for mixing tracer ffuctuations injected via steady bulk sources and sinks, where 
the dependence of the mixing measures on the the Peclet number depends explicitly on the 
length scales in the sources, the sinks, and those chosen for the mixing measure [3-5] . These 
questions remain for future investigations. 

Acknowledgments 

The authors are grateful to W. R. Young for helpful discussions, and to the hospitality of 
the Institute for Mathematics and its Applications (supported by NSF). CRD was supported 
by NSF under grants PHY-0555324 and PHY-0855335, and J-LT under grant DMS-0806821. 



[1] L. N. Howard, 'Heat transport by turbulent convection,' J. Fluid Mech. 17, 405 (1963). 

[2] J.-L. Thiffeault, C. R. Doering, and J. D. Gibbon, 'A bound on mixing efficiency for the 

advection-diffusion equation,' J. Fluid Mech. 521, 105-114 (2004). 
[3] C. R. Doering and J.-L. Thiffeault, 'Multiscale mixing efficiencies for steady sources,' Phys. 

Rev. E 74 (2), 025301(R) (2006). 
[4] T. A. Shaw, J.-L. Thiffeault, and C. R. Doering, 'Stirring up trouble: Multi-scale mixing 

measures for steady scalar sources,' Physica D 231 (2), 143-164 (2007). 



15 



[5] T. Okabe, B. Eckhardt, J.-L. Thiffcault, and C. R. Doering, 'Mixing effectiveness depends on 

the source-sink structure: Simulation results,' Journal of Statistical Mechanics: Theory and 

Experiment 2008 (7), P07018 (2008). 
[6] L. O Naraigh and J.-L. Thiffeault, 'Bounds on the mixing enhancement for a stirred binary 

fluid,' Physica D 237 (21), 2673-2684 (2008). 
[7] J.-L. Thiffcault and G. A. Pavliotis, 'Optimizing the source distribution in fluid mixing,' 

Physica D 237 (7), 918-929 (2008). 
[8] S. Plasting and W. R. Young, 'A bound on scalar variance for the advection-diffusion equa- 
tion,' J. Fluid Mech. 552, 289-298 (2006). 
[9] D. A. Birch, Y.-K. Tsang, and W. R. Young, 'Bounding biomass in the Fisher equation,' Phys. 

Rev. E 75, 066304 (2007). 
[10] S. R. Keating, P. R. Kramer, and K. S. Smith, 'Homogenization and mixing measures for a 

replenishing passive scalar field,' (2009), preprint. 
[11] M. R. Turner, J. Thuburn, and A. D. Gilbert, 'The influence of periodic islands in the flow 

on a scalar tracer in the presence of a steady source,' Phys. Fluids 21, 067103 (2009). 
[12] E. Balkovsky and A. Fouxon, 'Universal long-time properties of Lagrangian statistics in the 

Batchelor regime and their application to the passive scalar problem,' Phys. Rev. E 60 (4), 

4164-4174 (1999). 

[13] E. Gouillart, O. Dauchot, J.-L. Thiffeault, and S. Roux, 'Open-flow mixing: Experimental 

evidence for strange eigenmodes,' Phys. Fluids 21 (2), 022603 (2009). 
[14] E. Gouillart, O. Dauchot, and J.-L. Thiffeault, 'Measures of mixing quality in open flows with 

chaotic advection,' (2010), in preparation. 
[15] P. V. Danckwcrts, 'Continuous flow systems - distribution of residence times,' Chem. Eng. 
Set. 2, 1-13 (1953). 

[16] H. M. Hulburt, 'Chemical processes in continuous-flow systems,' Ind. Eng. Chem. 36 (11), 
1012-1017 (1944). 

[17] J. R. A. Pearson, 'A note on the "Danckwerts" boundary conditions for continuous flow 
reactors,' Chem,. Eng. Sci. 10, 281-284 (1959). 

[18] E. B. Nauman, 'Residence time distributions in systems governed by the dispersion equation,' 
Chem. Eng. Sci. 36 (6), 957-966 (1981). 

[19] A. Hisaka and Y. Sugiyama, 'Notes on the inverse Gaussian distribution and choice of bound- 
ary conditions for the dispersion model in the analysis of local pharmacokinetics,' J. Pharm. 
Sci. 88 (12), 1362-1365 (1999). 

[20] H. Aref, 'Stirring by chaotic advection,' J. Fluid Mech. 143, 1-21 (1984). 



